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The  main  drawback  of  Ni/YSZ  anode  supports  for  solid  oxide  fuel  cell  application  is  their  low  tolerance  to 
reducing  and  oxidizing  (RedOx)  atmosphere  changes,  owing  to  the  Ni/NiO  volume  variation.  This  work 
describes  a  structured  approach  based  on  design  of  experiments  for  optimizing  the  microstructure  for 
RedOx  stability  enhancement.  A  full  factorial  hypercube  design  and  the  response  surface  methodology 
are  applied  with  the  variables  and  their  variation  range  defined  as:  (1)  NiO  proportion  (40-60  wt%  of  the 
ceramic  powders),  (2)  pore-former  proportion  (0-30  wt%  corresponding  to  0-64vol.%),  (3)  NiO  particle 
size  (0.5-8  p,m)  and  (4)  8YSZ  particle  size  (0.6-9  pirn). 

To  obtain  quadratic  response  models,  25  different  compositions  were  prepared  forming  a  central  com¬ 
posite  design.  The  measured  responses  are  (i)  shrinkage  during  firing,  (ii)  surface  quality,  (iii)  as-sintered 
porosity,  (iv)  electrical  conductivity  after  reduction  and  (v)  expansion  after  re-oxidation.  This  approach 
quantifies  the  effect  of  all  factors  and  their  interactions.  From  the  quadratic  models,  optimal  compositions 
for  high  surface  quality,  electrical  conductivity  (>500  S  cm-1  at  room  temperature)  and  RedOx  expansion 
(<0.2%  upon  re-oxidation)  are  defined.  Results  show  that  expansion  after  re-oxidation  is  directly  influ¬ 
enced  by  the  sample  porosity  whereas,  surprisingly,  the  NiO  content,  varied  between  40  and  60  wt%,  does 
not  show  any  impact  on  this  response. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

As  energy  demand  is  increasing,  more  efficient  systems  and  sus¬ 
tainable  technology  are  required.  Solid  oxide  fuel  cells  (SOFCs)  can 
reach  high  efficiency  even  for  small  systems  and  are  able  to  use  a 
large  variety  of  fuels  (hydrogen,  natural  gas,  biogas,  etc.)  [  1  ].  Used  as 
electrolyzer,  they  can  store  excess  energy  production  from  renew¬ 
able  intermittent  sources  such  as  photovoltaics,  hydro-,  wind-  and 
tide-turbines. 

State-of-the-art  intermediate  temperature  (650-850  °C)  SOFC 
technology  is  based  on  a  thin  yttria-stabilized  zirconia  (YSZ) 
electrolyte  supported  by  a  nickel-YSZ  cermet  anode.  The  main 
drawback  of  this  design  is  its  sensitivity  to  nickel  re-oxidation 
(known  as  “RedOx”  cycling).  Oxidation  of  the  nickel  is  accompa¬ 
nied  by  a  volume  increase  of  71%  [2].  If  the  anode  support  and/or 
the  anode  active  functional  layer  (AFL)  expand  more  than  0.2%,  the 
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induced  tensile  stresses  will  provoke  cracks  in  the  thin  electrolyte 
[2-4].  The  re-oxidation  can  occur  (i)  at  high  current  density  or  high 
fuel  utilization,  (ii)  near  compressive  seals  [5],  (iii)  due  to  fuel  star¬ 
vation  and  (iv)  due  to  fuel  saving  during  startup  and  cool-down  of 
the  system.  The  RedOx  stability  of  the  cells  is  especially  important 
for  small  systems  (1-5  kW)  where  a  safety/emergency  gas  add-on 
is  not  economically  feasible. 

Changes  in  the  anode  microstructure  have  an  influence  on  the 
RedOx  behavior.  Waldbillig  et  al.  observed  that  the  expansion  upon 
RedOx  cycling  decreased  for  a  coarse  microstructure  compared  to 
a  fine  one  whereas  the  NiO  content  did  not  show  a  major  influ¬ 
ence  [6].  The  coarse  microstructure  enhances  the  RedOx  stability. 
In  a  similar  way,  Pihlatie  et  al.  showed  that  high  porosity  samples 
lower  the  expansion  during  RedOx  cycles  [7].  The  use  of  a  cermet 
containing  3  mol%  YSZ  showed  better  RedOx  stability  due  to  the 
higher  strength  of  3YSZ  compared  to  8YSZ  [3]. 

Microstructure  optimization  of  anodes  has  been  looked  at  since 
a  while  in  the  SOFC  community  [8].  These  optimizations  were 
applied  initially  to  increase  performance  by  enhancement  of  the 
electrical  conductivity,  the  porosity  and  triple  phase  boundary  den¬ 
sity  [9,10].  Later  the  focus  was  on  durability,  like  nickel  phase 
coarsening  problems  [11,12]. 
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Center  point  Factorial  point  Star  point 

Fig.  1.  Central  composite  design  for  three  factors  (/<  =  3). 


A  statistical  approach  for  optimization  of  a  complex  system  can 
reduce  the  number  of  experiments  compared  to  the  traditional  trial 
and  error  method.  This  approach  is  known  as  design  of  experiment 
(DoE)  and  response  surface  methodology  (RSM),  from  the  work  of 
R.A.  Fisher  in  the  first  quarter  of  the  last  century  [13,1 4].  This  math¬ 
ematical  method  gives  a  model  that  links  the  processing  conditions 
or  the  composition  of  a  product  (called  factors)  with  the  properties 
of  this  product  (measured  responses).  From  this  model,  optimal 
processes  and  mixtures  can  be  obtained  using  a  reduced  number 
of  experiments  [15,16].  The  approach  was  already  applied  to  SOFC 
anode  fabrication  for  plasma  spray  [17]  and  uniaxial  compaction 
methods  [18].  The  procedure  can  also  be  used  to  assess  the  main 
factors  responsible  for  long  term  degradation  [19]. 

The  goal  of  this  study  is  to  find  an  optimal  anode  support 
microstructure  in  particular  with  respect  to  re-oxidation  resis¬ 
tance.  The  desired  responses  of  the  support  are  (i)  an  expansion 
lower  than  0.2%  after  a  RedOx  cycle,  (ii)  sufficient  electrical  conduc¬ 
tivity  in  reduced  state  (>500  S  cm-1  at  room  temperature)  and  (iii) 
good  sinterability  (no  cracks  or  surface  defects  after  sintering).  The 
composition  of  the  prepared  samples  was  varied  using  different  fac¬ 
tors:  (1 )  the  NiO  content  in  the  ceramic  phase,  (2)  the  pore-former 
content,  (3)  the  NiO  particle  size  and  (4)  the  YSZ  particle  size. 

2.  Design  of  experiment 


a  hypercube  for  k>  3,  where  k  is  the  number  of  factors)  and  the 
levels  of  the  factors  are  defined  as  -1,  0,  +1  (for  the  minimal,  the 
intermediate  and  the  maximal  values  of  the  factor). 

Response  surface  methodology  requires  a  model  for  the  true 
response  surface  approximation.  If  no  physical  model  is  known, 
multiple  regressions  can  be  used  to  build  an  empirical  model.  The 
commonly  used  technique  is  the  least  squares  regression  model: 

y  =  Xa  +  e  (1) 


where  y  is  the  vector  of  the  n0  observations,  a  is  the  vector  of  the 
rip  regression  coefficients  parameters,  X  is  the  n0  x  np  model  matrix 
built  from  the  design  of  experiment  and  e  is  the  vector  of  random 
errors.  The  model  matrix  (based  on  -1  and  +1 )  will  be  converted  into 
the  experimental  matrix  when  the  range  of  the  factors  is  defined  (see 
next  section  and  Table  1).  The  unbiased  least  squares  estimator  of 
a  is  given  by 

a  =  (X'X)_1X'a  (2) 

where  X'  is  the  transposed  model  matrix,  (X'X)-1  is  the  dispersion 
matrix.  The  built  model  is  then  given  by 


4  4  4 

Y(x)  =  a0  +  £**+££  QgXjX;  +  e(0,  cr2)  (3) 

1=1  j=  1  1=1 


with  Y(x)  an  estimation  of  the  measured  response^,  a0  the  constant 
effect,  a*  the  main  half-effects  of  factors  x2,  ^  the  first  order 

interaction  half-effects,  a,-,-  the  quadratic  half-effects  and  £(0 ,<r2) 
the  residuals  that  should  follow  a  Gaussian  distribution  centered 
on  zero  (E(s)  =  0)  with  a  variance  of  o2  ( Var(e)  =  o2 ). 

The  comparison  between  different  models  of  the  same  response 
(i.e.  linear  model  versus  quadratic  model)  can  be  made  using  the 
adjusted  coefficient  of  linear  regression,  R2dj ,  as  the  standard  coeffi¬ 
cient  of  linear  regression,  R2,  will  always  improve  when  parameters 
are  added  to  the  model: 

R2  =  1  -  ^  (4) 


Kis-l-^O-K2) 


n0  -  n» 


(5) 


where  SSj  is  the  sum  of  squares  of  the  observations,  and  SSe  is  the 
sum  of  squares  of  the  errors  or  residuals  (the  residual  i  is  the  dif¬ 
ference  between  observation,  y2  and  the  fitted  value,  Y2),  n0  is  the 
number  of  observations  and  np  is  the  number  of  parameters  used 
in  the  model. 


2. 1.  Factorial  design  and  response  surface 


2.2.  Experimental  factors  and  measured  response 


The  idea  of  factorial  design  is  to  observe  the  effects  of  the  factors 
and  their  interaction  on  the  measured  responses  with  relatively 
fewer  experiments.  Factors  can  be  quantitative  (e.g.  concentration) 
or  qualitative  (e.g.  different  catalysts  or  different  powder  types). 
The  simplest  design  uses  only  two  levels  of  each  factor,  basically 
the  minimum  and  the  maximum,  which  are  noted  in  coded  value 
-1  and  +1,  respectively,  for  simplification.  For  full  factorial  design 
all  the  possible  combinations  are  executed,  but  when  the  factors 
are  numerous,  fractional  factorial  design  can  be  used  [15,16]. 

The  central  composite  design  (CCD)  is  widely  used  for  response 
surface  methodology  (RSM)  for  a  second-order  model  [15,16].  This 
design  involves  two  overlapping  designs:  a  two  level  factorial 
design  represented  by  the  open  symbols  in  Fig.  1  and  an  axial  or 
star  design  represented  by  the  stars  in  Fig.  1.  The  factorial  points 
represent  the  variance-optimal  design  for  a  first  order  plus  two- 
factor  interaction  model.  The  center  point  and  star  points  allow 
efficient  estimation  of  the  quadratic  terms  (see  [16],  p.  297).  With 
a ,  the  distance  of  the  star  points  from  the  center  point,  equal  to 
1,  the  star  points  are  placed  in  the  center  of  the  cubic  faces  (or  of 


The  factors  are  the  studied  experimental  variables  that  could 
influence  the  sample  properties  (also  noted  as  measured  responses 
of  the  sample).  The  experimental  factors  in  this  study  are:  (1)  the 
NiO  weight  fraction  in  the  ceramic  phase,  (2)  the  pore-former  con¬ 
tent  added  to  the  ceramic  powders,  (3)  the  NiO  particle  size  and  (4) 
the  8YSZ  particle  size.  All  the  natural  and  coded  values  are  sum¬ 
marized  in  Table  1.  The  natural  value  is  the  chosen  interval  for  the 
factor  (e.g.  0, 1 5  and  30  wt%  of  pore-former)  and  the  corresponding 
coded  (or  standardized)  variables  change  from  -1  to  +1  on  the  full 
interval.  The  use  of  the  coded  values  is  for  the  ease  of  the  matrix 
calculus. 


2.2  A.  Weight  fraction  of  NiO 

The  weight  fraction  of  NiO  is  calculated  only  for  the  ceramic 
phase  (i.e.  excluding  pore-former).  The  natural  variable  for  NiO 
weight  fraction  (u1Ni0)  is  obtained  as  follows: 


ul,NiO 


mNiO 

^Nio  +  msz 


x  100 


(6) 
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Table  1 

Variables  or  factors  studied:  from  natural  to  coded  values. 


Factors\coded  value  x, 

-1 

0 

1 

U\ 

NiO /YSZ  fraction 

40/60  wt% 

50/50  wt% 

60/40  wt% 

U2 

Pore-former  fraction 

0wt% 

15  wt% 

30  wt% 

u3 

NiO  particle  size 

Fine 

50  wt%  fine/50  wt%  coarse 

Coarse 

u4 

8YSZ  particle  size 

Fine 

50  wt%  fine/50  wt%  coarse 

Coarse 

with  40<u1Nio  <60,  mNio  and  mYsz  the  mass  of  NiO  and  YSZ, 
respectively. 

The  standardized  variable  for  NiO  weight  fraction  {x^  Ni0)  is 

v  ul,NiO  -  ul,Nio(0)  .  . 

xl,NiO  =  - ^ -  V') 

with  ii!  Nio(0)  the  center  of  the  natural  interval  and  A u  half  of 
the  natural  interval.  Then  the  interval  for  the  natural  variable 
^i,Nio  =  [40,60]  corresponds  to  the  interval  for  the  standardized 
variable  Ni0  =  [-1,1  ].  This  is  done  for  all  the  measured  variables. 
The  natural  variable  for  YSZ  weight  fraction  (ui  ysz)  is 

Ml, YSZ  =  100  -  Ui  NiO  (8) 


15  min  of  applying  the  ultrasound  horn.  The  refractive  indexes 
used  for  the  powders  are  2.37,  2.10  and  1.52  for  NiO,  YSZ  and 
pore-former,  respectively  and  the  densities  are  1.55  gem-3  for  the 
pore-former  and  summarized  in  Table  B1  for  NiO  and  YSZ  [20,21  ]. 

Specific  surface  area  measurements  were  carried  out  with  a 
Micrometries  Gemini®  2376  instrument  using  a  classical  5  point- 
measurement.  Prior  to  the  measurement,  the  samples  were  heated 

1  h  at  200  °C  under  a  nitrogen  flux  to  remove  any  adsorbed  water. 

Thermogravimetric  analysis  was  performed  with  a  TGA 
SDTA851  instrument  from  Mettler  Toledo®  at  a  ramp  rate  of 

2  °C  min-1  from  room  temperature  up  to  800  °C  under  air  flow  in  a 
150  pul  alumina  crucible. 


2.2.2.  Pore- f ormer  (PF) fraction 

The  weight  proportion  of  pore-former  (u2)  is  added  to  the  total 
weight  of  ceramics  (NiO  +  YSZ): 


u2 


™ Pore-former 

mN  iO  +  mYSZ  +  ™ Pore-former 


X  100 


withO  <u2  <30 


(9) 


Using  equations  in  Appendix  B  and  densities  in  Table  B1  and 
1.55 gem-3  for  the  pore-former  (PF),  15wt%  PF  corresponds  to 
41  -42  vol.%  PF,  depending  of  the  weight  fraction  of  NiO  and  30  wt% 
PF  corresponds  to  63-64  vol.%  PF. 


2.2.3.  NiO  and  8YSZ  particle  size 

Two  different  types  of  nickel  oxide  and  8YSZ  were  used  here.  The 
powder  characterization  will  be  described  later.  The  standardized 
variable  for  NiO  size  is  given  by 

X3=2 - - 1  (10) 

™fine  +  ^coarse 

The  standardized  variable  for  the  8YSZ  particle  size  is  calculated 
in  the  same  way. 

The  measured  responses  are  the  following: 

•  Shrinkage  after  firing  4  h  at  1400  °C  in  air. 

•  Qualitative  observation  of  the  sample  surface  after  sintering. 

•  Porosity  in  as-sintered  and  reduced  state. 

•  Conductivity  at  room  temperature  after  reduction  for  24  h  at 
800  ±  5  °C  in  a  flow  of  8%  H2  in  N2. 

•  Expansion  at  room  temperature  after  re-oxidation  during  24  h  at 
845  ±  5  °C  in  air. 


3.  Experimental 

3.1.  Raw  powder  characterization 

Different  powders  were  used  for  this  study,  two  NiO  and  two 
8YSZ  (8mol%  Y203-Zr02)  powders,  and  an  organic  pore-former. 
Particle  size  measurements  were  carried  out  by  laser  diffraction 
and  specific  surface  area  determination.  The  pore-former  was  also 
studied  by  thermogravimetric  analysis  (TGA). 

Particle  size  distributions  were  obtained  by  laser  diffraction 
using  a  Mastersizer  (Malvern  Inst.  (UK)),  with  Mastersize  S  2.15 
software.  The  analysis  model  is  polydisperse  and  the  refractive 
index  of  the  solvent  is  1.33.  The  analysis  was  done  in  an  aque¬ 
ous  suspension  with  0.02%  PAA  (poly  acrylic  acid)  dispersant  after 


3.2.  Sample  preparation  and  microstructure  comparison 

Planar  SOFC  anode  supports  are  usually  fabricated  by  tape  cast¬ 
ing,  a  well-known  technique  easy  to  upscale  [22,23].  The  technique 
uses  slurries  dispersed  by  a  wet  ball-milling  process  with  an  opti¬ 
mized  composition  of  solvent,  dispersant,  binder,  plasticizer  and 
defoamer.  As  this  study  is  focused  on  the  powder  size,  pore-former 
and  the  phase  proportions  of  the  samples,  and  not  on  the  opti¬ 
mization  of  the  wet  slurry,  a  dry  ball-milling  process  was  used, 
with  Turbula®  equipment  followed  by  an  isostatic  compaction  in 
cylinder  shape. 

The  duration  of  dry  ball-milling,  with  0.5  wt%  of  stearic  acid 
as  dispersant,  was  optimized  to  obtain  the  same  microstructure 
than  after  standard  wet  ball-milling  during  24  h.  The  microstruc¬ 
ture  after  2,  4,  6,  8  and  10  h  of  dry  ball-milling  was  compared  to 
standard  samples  (obtained  from  24  h  wet  ball-milling)  by  using 
a  scanning  electron  microscope  (FEI  Xlf-30  SFEG  Sirion,  at  lOkV, 
with  backscattered  electrons  (BSE)  detector,  spot  size  4,  work¬ 
ing  distance  5  mm  and  a  magnification  of  1500x)  after  polishing 
the  NiO-YSZ  samples  down  to  0.5  p,m  with  diamond  lapping  film 
and  water  as  lubricant  and  cooling  medium.  From  these  obser¬ 
vations,  8  h  of  dry  milling  in  a  Turbula®  mill  appeared  to  yield 
a  microstructure  closest  to  the  standard  one.  For  shorter  times, 
agglomerates  remained,  whereas  for  longer  times,  the  microstruc¬ 
ture  became  finer  due  to  further  particle  size  reduction  of  the 
powders.  Therefore,  the  dry  ball-milling  process  for  the  tested 
samples  was  fixed  at  8  h  in  the  Turbula®  mill  with  0.5  wt%  of 
stearic  acid  (relative  to  total  powder  mass)  and  1  more  hour 
with  1  wt%  of  Duramax®  B1000  acrylic  binder  (relative  to  total 
mass). 

To  obtain  homogeneous  microstructures  in  a  compacted  sam¬ 
ple,  the  mixed  powder  was  isostatically  pressed  at  200  MPa  for 
1  min  (ramp  of  1  kN  s-1  and  plateau  at  400  kN)  in  a  silicone  cylinder 
of  10  mm  diameter  and  40  mm  length.  Sample  faces  were  polished 
parallel  with  600  mesh  SiC  paper,  for  their  length  measurement 
done  with  a  micrometer. 

Sintering  was  done  for  4h  in  air  at  1400  °C,  with  a  1  °Cmin-1 
ramp  from  room  temperature  to  800  °C  and  a  3  °C  min-1  ramp  from 
800  to  1400  °C.  Cooling  occurred  with  a  1  °C  min-1  ramp  from  1400 
to  1 100  °C  and  then  naturally  down  to  room  temperature. 

Reduction  was  carried  out  in  a  tubular  furnace.  The  procedure 
was  to  first  heat  at  10°Cmin-1  up  to  800°C±5°C  under  an  air 
flux  of  1  lmin-1.  After  1.5  h  of  temperature  stabilization,  nitrogen 
was  introduced  for  20  min  at  1 .8 1  min-1  for  a  purge.  Then  8%  FI2  in 
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Measurements  taken  on  the  cylindrical  samples  at  the  different  stages. 
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Before  sintering 

After  sintering 

After  reduction 

After  re-oxidation 

Mass  [g]  ±0.001  g 

V 

V 

V 

V 

Length  [mm]  ±0.001  mm 

V 

V 

V 

V 

Diameter  [mm]  (two  measurements  at  the  extremities  and  two  at  the  center)  ±  0.01  mm 
Conductivity  (S  cm-1 ) 

V 

V 

Table  3 

Powder  characterizations. 

Powders 

Light  diffraction,  dv 50  (pm)  (span) 

Specific  surface  area,  Sbet  (m2  g-1 ) 

deET  (pm) 

Factor  of  agglomeration  d^o/deET 

NiO  fine 

0.52  (2.88) 

4.29 

0.209 

2.5 

NiO  coarse 

8.37(2.25) 

2.79 

0.322 

26.0 

8YSZ  fine 

0.57  (3.45) 

10.98 

0.090 

6.3 

8YSZ  coarse 

9.08(1.07) 

0.66 

1.502 

6.0 

Pore-former 

8.56(1.61) 

0.78 

4.962 

1.7 

N2  (form  gas)  was  flown  at  1-1.51  min-1  for  24  h,  and  finally  the 
furnace  cooled  down  at  10°Cmin-1  under  form  gas. 

A  similar  procedure  was  applied  for  oxidation,  switching  gas 
appropriately.  The  plateau  was  fixed  at  845  ±  5  °C.  This  error  was 
due  to  the  temperature  distribution  inside  the  tubular  furnace. 
Temperature  increased  when  starting  the  air  flow,  due  to  the  strong 
exothermicity  of  nickel  oxidation.  During  the  experimental  proce¬ 
dures,  sample  fabrication  was  reproduced  from  1  to  4  times. 


3.3.  Sample  measurements 

Measurements  taken  on  the  cylindrical  samples  at  the  different 
stages  are  listed  in  Table  2.  As  no  variation  of  length  was  measured 
between  as-sintered  and  reduced  state,  the  shrinkage  of  the  sample 
during  reduction  is  neglected. 

Porosity  proportion  (ppore)  of  the  sample  was  extracted  directly 
from  the  measured  mass  (msamp/e)  and  volume  of  the  sample 

sample )• 


„  _  VPore 

Ppore  -  w 

v  sample 


1  tn sample 
Pth  y sample 


(ii) 


with  pth  density  of  the  composite  materials  (known  from  the  com¬ 
position). 

Conductivity  (or)  in  S  cm-1  of  a  sample  is  obtained  by  measuring 
the  potential  drop  (AV)  in  V  between  two  contacts  separated  by  a 
known  distance  (d)  in  cm: 


a[Scm"11  =  ^i7  (12) 

with  r  the  radius  in  cm  of  the  sample  and  I  a  constant  current  in  A 
passed  through  the  sample  (4-point  method). 


4.  Results  and  discussion 

4.1.  Powder  characterizations 

A  summary  of  the  powder  characterizations  is  given  in  Table  3 
The  dv50  is  the  equivalent  volume  spherical  diameter  at  which  half 
of  the  particles  volume  is  bigger  and  half  is  smaller  than  this  value 
[24].  The  span  gives  a  measure  of  the  particle  size  distribution  and 
is  given  by 

Cn,n  dvgo  -  dyio 

Span  = - - -  (13) 

50 

where  dv io  and  dv 90  are  particle  sizes  where  10  and  90vol.%  are 
smaller  than  this  value,  respectively. 


From  the  specific  surface  measurement  (SBET),  a  monodisperse 
spherical  diameter  (c/bet)  can  be  calculated: 


dBET[fxm]  = 


_ 6 _ 

SbetIhi2  g-1]  p[gcm-3] 


(14) 


where  p  is  the  powder  density. 

The  light  scattering  results  show  that  the  fine  powders  have  a 
diameter  of  about  half  a  micron  and  that  the  coarse  powders  and  the 
pore-former  have  a  diameter  around  9  pm.  But  if  we  compare  the 
specific  surface  area  measurements,  the  coarse  NiO  shows  agglom¬ 
erations  that  are  not  broken  down  during  the  ultrasound  treatment 
before  light  scattering  measurements.  These  agglomerates  are  also 
observed  in  Fig.  2b.  The  energy  of  the  ball-milling  treatment  is 
higher  than  the  one  from  the  ultrasound  treatment,  so  agglom¬ 
erates  can  break  during  the  homogenizing  with  ball-milling.  Fig.  2 
shows  scanning  electron  micrographs  of  the  different  powders. 

Fig.  3  shows  the  TGA  measurement  of  the  pore-former.  The  first 
decrease  to  90  wt%  is  due  to  the  evaporation  of  the  water  included 
in  the  pore-former.  It  is  followed  by  a  stronger  decrease  around 
300  °C  corresponding  to  the  release  of  the  chemically  bonded  water. 
The  ultimate  mass  variation  from  350  to  600 °C  comes  from  the 
decomposition  of  the  organic  phase. 


4.2.  Raw  results 

As  some  sample  compositions  were  replicated  the  total  number 
of  samples  is  46  from  25  different  compositions  (at  least  one  sample 
of  each  composition  was  tested).  The  central  point  was  repeated 
three  times.  The  full  list  of  the  25  different  compositions  and  the 
raw  results  are  presented  in  supplementary  on-line  data. 

4.2.1.  Expansion  after  re -oxidation 

The  results  show  a  certain  reproducibility  depending  on  the 
composition  of  the  sample.  The  main  goal  of  this  study  being  to 
obtain  samples  with  an  expansion  after  a  RedOx  cycle  lower  than 
0.2%,  one  sees  that  samples  from  compositions  1,  3,  9,  11  and  20 
fulfill  this  condition.  It  is  interesting  to  note  that  all  these  samples 
contain  only  coarse  YSZ  powder. 

The  sample  surface  can  be  without  any  apparent  defects  (surface 
quality,  SQ  =  3),  show  some  holes  and  small  cracks  (SQ  =  2)  or  larger 
cracks  (SQ=  1)  (see  Fig.  4a-c).  The  low  expansion  samples  are  all 
poor  in  surface  quality  as  shown  by  Fig.  4d.  Only  the  composition 
3  could  give  an  intermediate  SQ  value.  A  low  SQ  value  will  not  be 
appropriate  from  a  fabrication  point  of  view. 

Expansion  after  re-oxidation  is  strongly  influenced  by  the  poros¬ 
ity  of  the  samples  as  shown  in  Fig.  5.  With  as-sintered  porosity 
lower  than  35  vol.%,  it  is  not  possible  to  achieve  the  0.2%  expansion 
limit.  Flowever,  also  several  samples  with  porosity  above  45  vol.% 
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Fig.  2.  Scanning  electron  micrographs  of  the  powders:  (a)  NiO  fine,  (b)  NiO  coarse,  (c)  8YSZ  fine,  (d)  8YSZ  coarse  and  (e)  pore-former  (marker  length  =  5  p,m  for  all  micrographs). 


showed  a  slightly  too  high  expansion.  The  relation  between  expan¬ 
sion  after  re-oxidation  (e)  and  the  as-sintered  porosity  (ppore)  could 
be  approximated  with  an  exponential  or  power  law: 


-exp  —  8.825  exp(— 0.081  -Ppore) 

(15) 

'Pow  =  697.34(ppore) 

(16) 

The  R 2  coefficient  was  slightly  better  for  the  exponential  fit 
(^exp  =  0.725)  compared  to  the  one  from  the  power  law  (flp0W  = 
0.693).  The  0.2%  RedOx  expansion  limit  corresponded  to  47  and 
54%  porosity  proportion  for  exponential  and  power  law  approx¬ 
imations,  respectively.  In  other  words,  a  NiO-YSZ  anode  support 
needs  a  porosity  after  sintering  around  50%  to  be  RedOx  stable. 

Pihlatie  et  al.  showed  that  changing  the  porosity  of  a  sample 
from  9  to  33%  by  reducing  the  temperature  of  firing  by  20  K  and 
introducing  pre-calcined  NiO  powder,  the  cumulative  RedOx  strain 
was  decreased  by  a  factor  of  20  [7].  If  we  introduce  these  porosity 
values  into  Eqs.  (15)  and  (16),  the  resulting  factor  decrease  in  linear 
expansion  is  calculated  to  be  equal  to  7  and  14,  respectively.  The 
trend  in  the  present  study  is  thus  quite  similar  compared  to  the 
findings  of  Pihlatie  et  al.  [7]. 

4.2.2.  Porosity  after  reduction 

For  convenience,  the  composition  of  a  cermet  is  prepared  in 
weight  percentage  (wt%),  w,.  Conductivity  calculation  is  related  to 
the  volume  percentage  (vol.%)  of  the  conducting  phase  in  the  solid 
(vNi)  as  well  as  the  porosity  fraction  (pp0re).  The  weight  and  the  vol¬ 
ume  fractions  in  the  cermet  change  during  the  chemical  reaction 
of  reduction  or  oxidation  due  to  the  variation  of  molar  mass  and 
molar  volume.  The  molar  percentage  (mol%),  n2-,  of  Ni  as  NiO  or  Ni 
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Fig.  6.  Oxidized  porosity  (as  sintered)  versus  reduced  porosity  (24  h  at  800  °C  in 
8mol%  H2  in  N2)  measurements  (diamonds:  60wt%  NiO,  triangles:  50wt%  NiO, 
crosses:  40wt%  NiO)  compared  to  Eq.  (17)  (lines). 


Fig.  4.  Example  of  the  different  surface  qualities  (SQ):  (a)  SQ  =  1  (composition  3), 
(b)  SQ  =  2  (composition  21),  (c)  SQ  =  3  (composition  16),  and  (d)  expansion  after 
re-oxidation  versus  SQ.(the  dashed  line  is  the  RedOx  limit).  Marker  length  =  1  cm. 


stays  constant  in  the  cermet.  The  fraction  of  the  nickel  or  nickel 
oxide  phase  in  the  composite  is  given  in  Table  4  for  the  different 
studied  samples.  The  equations  and  the  values  used  to  construct 
the  table  are  given  in  Appendix  A. 

If  the  volume  of  the  sample  is  constant  during  reduction,  the 
porosity  variation  can  be  calculated  by  combining  Eqs.  (B4),  (Bll) 
and  (B12).  This  gives  the  following  relation  between  the  initial 
porosity  after  sintering  (oxidized  state),  pp0re,ox  and  the  reduced 
sample  porosity,  ppore,red • 


PPore,red  —  f  +  (1  —  £)Ppore,ox 


(17) 


with  £  = 


1-(WW 

1  +  ((1/w)  -  l)(PNio/PYSz) 


where  w  is  the  weight  fraction  of  NiO,  VNi  and  VNi0  are  the  molar 
volume  of  nickel  and  its  oxide,  respectively,  pNi0  and  pY sz  are 
the  respective  densities.  Radovic  and  Lara-Curzio  had  a  similar 
approach  using  a  different  equation  [25]. 

The  measured  porosity  before  and  after  reduction  for  24  h  at 
800  °C  in  8%  H2  in  N2  are  compared  to  equation  17  in  Fig.  6.  The 
model  is  accurate  for  the  40wt%  NiO  measurements,  but  slightly 
overestimates  the  50  and  60wt%  NiO  values.  This  could  be  due  to 
the  fact  that  reduction  is  not  completed  after  24  h  for  a  cermet  with 
higher  nickel  content.  As  the  samples  are  relatively  thick  (diameter 
between  5  and  8  mm  depending  on  the  sintering  shrinkage),  the 
reduction  will  take  a  certain  time,  as  the  gas  has  to  diffuse  in  the 
entire  volume.  The  higher  the  nickel  oxide  content,  the  longer  the 
time  needed  for  complete  reduction.  The  volume  fraction  of  NiO 
varies  from  37  to  57  vol.%  NiO,  the  nickel  oxide  can  agglomerate  and 
forms  coarse  particles  after  sintering  for  high  NiO  fraction  samples. 


4.2.3.  Electrical  conductivity 

Conductivity  of  the  sample  will  depend  on  the  proportion  of  the 
different  conducting  phases  as  well  as  on  the  porosity  fraction.  The 
basic  conductivity  approximation  is  obtained  from  the  mean  value 
of  the  conductivities  as  a  function  of  the  phase  fractions  (model  1 ): 

<JcomP,  1  =  PNi^Ni  +  PYSZ^YSZ  +  Ppore<7pore  (18) 

where  are  the  electrical  conductivities  of  the  composite  and  the 
different  phases  and  p*  are  the  volume  fractions  of  each  phase  in 
the  sample  (pNi  +pY sz  +Ppore  =  1).  The  electrical  conductivity  of  Ni 
and  8YSZ  at  25  °C  is  equal  to  146,200  and  0.006  S  cm-1 ,  respectively 
[9,26].  As  the  conductivity  of  the  pore  is  assumed  equal  to  zero,  the 


Table  4 

Comparison  of  the  nickel  or  nickel  oxide  fraction  for  the  studied  sample  (w:  weight  fraction,  v:  volume  fraction  in  the  solid  phase,  n:  mole  fraction,  p,:  volume  fraction 
including  porosity,  and  £,  the  origin  of  Eq.  (17)). 


Weight  fraction  oxidized  :wNi0 

40wt%  NiO 

50  wt%  NiO 

60wt%  NiO 

Molar  fraction:  nNi0  =  nNi 

0.533 

0.631 

0.720 

Volume  fraction  oxidized:  vNio 

0.371 

0.469 

0.570 

Volume  fraction  reduced:  vNi 

0.258 

0.343 

0.439 

Weight  fraction  reduced  :wNi0 

0.344 

0.440 

0.541 

£  (reduced  porosity  for  dense  oxidized  sample,  pp0re,ox  =  0) 

0.152 

0.193 

0.234 
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Fig.  7.  Linearization  of  model  3  (Eq.  (22)).  Measurements  in  the  dashed  circle  could 
be  due  to  bad  electrical  contacts  during  measurement. 


•  60wt%NiO  (6N) 

°  50wt%NiO  (5N) 

*  40wt%NiO  (4N) 

— Model  1  (6N) 

— Model  1  (4N) 

— Model2,  k=3  (6N) 

—  -Model2,  k=3  (4N) 
- Model3,  q=1.83  (6N) 

—  Model3,q=1.83  (5N) 

—  -Model3,  q=l  .83  (4N) 


Fig.  8.  Conductivity  versus  the  solid  fraction  (pso/id  =  ^  -  Ppore )  and  the  NiO  weight 
fraction  in  the  ceramic  mixture.  Measurements  in  the  dashed  circle  could  be  due  to 
bad  electrical  contact  during  measurement  (same  measurements  than  in  Fig.  7). 


last  equation  can  be  rewritten  as 

G comp,  1  =  (1  ~  Ppore)&eff  (19) 

With  crejj  =  UNi^Ni  +  WSZ^YSZ 

where  uNi  and  uYsz  are  the  Ni  and  the  YSZ  volume  fraction  in  the 
solid  phase  (vNi  +  uYsz  =  1  and  1  -ppore  =PNi/%i  =Pysz/vYsz)-  A  sim¬ 
ilar  model  uses  an  additional  power  law  parameter,  /<,  to  fit  the 
conductivity  [27]  (model  2)  to: 

&comp,2  —  (1  ~  Ppore)  &eff  (20) 

where  k  should  take  a  value  between  1.5  and  3. 

Here,  we  propose  a  model  including  a  k  value  varying  with  the 
conducting  phase  fraction  (vNi);  Eq.  (20)  is  changed  to  model  3: 

° comp, 3  —  (1  ~  Ppore)^/  Nl&eff  (21) 

The  q  value  can  be  extracted  by  linearization  of  Eq.  (21): 

ln(crC0mp?3)  —  \n(creff)  =  —  ln(l  —  pp0re )  (22) 

•'Ni 

Fig.  7  shows  the  plots  of  Eq.  (22)  with  our  experimental  values  for 
model  3.  With  this  model  the  different  NiO  fractions  are  grouped 
together.  The  observation  shows  that  four  values  lie  far  outside  of 
this  linear  approximation  (the  four  points  in  the  dashed  circle).  The 
hypothesis  is  that  these  values  were  due  to  bad  electrical  contacts. 
Without  these  measurements,  the  R2  coefficient  of  the  lineariza¬ 
tion  increases  from  0.9  to  0.97  and  the  q  value  decreases  from  1.98 
to  1.83.  The  latter  q  value  is  used  to  fit  the  experimental  measure¬ 
ments  in  Fig.  8  where  the  different  models  are  compared. 

The  simplest  model  1  shows  much  larger  values  than  the  mea¬ 
surements.  Model  2  with  the  exponent  /<  =  3  gives  a  reasonable 
upper  bound  for  the  samples  containing  60wt%  NiO  (xi=+l). 
The  conductivity  measurements  of  the  samples  containing  50  and 
40wt%  NiO  are  more  dispersed  and  no  value  lies  close  to  the 
upper  limit  given  by  model  2.  Values  for  two  samples  with  40  wt% 
NiO  were  equal  to  zero  and  not  shown  in  Fig.  8.  The  conductiv¬ 
ity  is  strongly  influenced  by  the  conducting  phase  fraction  (vNi), 
as  shown  by  the  mean  variation  with  the  NiO  content:  from  5250 
to  2680  and  11 30 Son-1  for  60,  50  and  40wt%  NiO,  respectively. 
Model  3,  proposed  here  (Eq.  (21 )),  is  incorporating  the  influence  of 
the  volume  fraction  of  conducting  phase  on  the  electrical  conduc¬ 
tivity.  A  few  points  lie  outside  of  the  model  3  approximation:  this 
can  be  due  again  to  a  contact  issue  in  the  measurement  method, 
since  it  appears  for  the  three  different  NiO  compositions.  But  it  can 
also  come  from  a  percolation  issue,  as  two  samples  could  not  be 
measured  at  all  (sample  10  and  16,  see  additional  on-line  support). 


The  measured  samples  contain  a  relatively  small  nickel  fraction 
if  the  porosity  volume  is  included,  pNi  =  vNi(l  -  Ppore ),  values  rang¬ 
ing  from  13.9  to  32.4vol.%  Ni.  The  13.9vol.%  Ni  sample  (number 
9)  showed  a  conductivity  of  1 70S  cm-1  while  the  non-conductive 
samples  (10  and  16)  had  a  pN j  of  16.1  and  18.1  vol.%  Ni,  both  orig¬ 
inally  containing  40wt%  NiO.  Costamagna  et  al.  showed  that  the 
percolation  threshold  in  a  two-phase  dense  packed  sphere  model 
depends  on  the  particle  size  ratio  (dNi/dYSz)  [28].  If  the  ratio  is  close 
to  unity,  the  volume  fraction  for  percolation  is  0.294.  If  the  ratio  is 
equal  to  0.3  then  the  percolation  should  occur  around  0.1 2.  Conduc¬ 
tion  should  be  possible  even  for  low  amounts  of  nickel  if  its  particles 
are  (much)  smaller  than  YSZ  or  smaller  than  the  porosity.  This 
model  is  neither  taking  into  account  the  high  amount  of  porosity 
nor  the  fact  that  the  real  particle  shapes  are  different  from  spheres. 
But  we  can  consider  that  porosity  is  a  “non-electronically  conduc¬ 
tive  medium”  similar  to  YSZ.  Then  sample  9  with  pNi  =  13.5  vol.% 
could  have  had  a  ratio  between  nickel  particles  and  YSZ  and  pore 
particles  smaller  than  1/3.  It  is  important  to  note  that  the  con¬ 
ductivity  measurements  are  done  after  24  h  under  forming  gas  at 
800  °C:  it  is  well  known  that  the  nickel  phase  can  sinter  under  nor¬ 
mal  SOFC  conditions  (high  temperature  and  high  vapor  pressure) 
[11,12].  After  nickel  coarsening,  the  particles  ratio  between  Ni  and 
YSZ  will  increase  and  the  percolation  of  the  Ni  phase  could  be  lost. 

4.3.  Surface  response  and  model  comparison 

Comparison  of  the  linear  and  quadratic  models  is  given  in  Table  5 
using  the  adjusted  coefficient  of  multiple  determination  R2dj  (Eq. 
(5)).  In  all  cases  the  quadratic  model  is  better,  especially  for  the 
conductivity  and  expansion  responses.  If  the  addition  of  parameters 
between  linear  and  quadratic  models  did  not  enhance  significantly 
the  fit,  the  value  of  R2dj  would  have  decreased.  The  model  used  for 
the  response  estimation  is  thus  in  all  cases  quadratic.  The  regression 
coefficients  of  the  different  responses  are  given  in  Table  6.  The  a0 
estimator  is  the  mean  value  of  the  full  set  of  measured  samples.  For 
example,  the  mean  value  for  expansion  is  1.39%,  a2  equal  to  -1.04 
means  the  expansion  will  decrease  by  -1 .04  when  the  pore-former 
weight  content  changes  from  u2  =  15-30%,  corresponding  to  x2  =  0 
and  1,  respectively.  When  the  interaction  effect  is  positive,  e.g.  ai4, 
the  response  value  will  increase  if  factors  1  and  4  are  both  positive 
or  both  negative.  Negative  quadratic  effects,  such  as  a33  and  a44, 
will  give  a  maximum  inside  the  tested  region  and  minimal  values 
at  the  border  as  shown  by  Fig.  9a.  The  result  of  high  values  of  a2 
and  a22  compared  to  low  values  of  a\,  an  and  ai2,  shows  that  the 
response  surface  for  the  expansion  model  depends  mainly  on  the 
pore-former  content  (x2)  in  Fig.  9b,  and  hardly  on  NiO  content. 
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Table  5 


Adjusted  coefficient  of  linear  regression,  R^dj  (Eq.  (6))  for  linear  and  quadratic  models  of  the  different  responses. 

R^d.  Expansion  after  oxidation  (%)  Conductivity  (S  cm-1 )  Surface  quality 

Porosity  (%) 

Shrinkage  (%) 

Linear  model  0.792  0.748  0.949 

Quadratic  model  0.861  0.879  0.962 

0.973 

0.989 

0.995 

0.996 

Table  6 

Least  square  estimates  of  the  quadratic  models  for  the  different  measured  responses. 

Estimate 

Expansion  (%) 

Conductivity  (S  cm-1 ) 

Surface  quality 

Porosity  (%) 

Shrinkage  (%) 

«o 

1.39 

2124 

2.29 

30.1 

-16.9 

Qi 

0.10 

2455 

-0.15 

1.8 

0.4 

a2 

-1.04 

-1364 

-0.53 

9.5 

-6.5 

a3 

0.30 

794 

-0.04 

-1.9 

-0.5 

a4 

-0.56 

330 

-0.32 

4.5 

3.6 

an 

-0.03 

-738 

-0.01 

2.0 

0.3 

an 

-0.04 

1349 

0.03 

-1.0 

-0.1 

au 

0.31 

-271 

0.17 

-0.2 

-0.4 

<223 

-0.27 

-1553 

-0.12 

1.9 

0.4 

c*24 

0.07 

-1034 

-0.31 

2.6 

0.4 

a34 

0.03 

120 

-0.05 

1.3 

-0.1 

an 

-0.26 

18 

-0.11 

2.1 

1.0 

a2  2 

1.53 

1090 

0.28 

-5.7 

-1.9 

033 

-0.44 

258 

0.48 

-3.0 

-0.4 

a4  4 

-0.76 

-211 

-0.61 

1.8 

0.9 

NiO  particles 
size 


(b) 


Expansion  % 


Ratio  NiO/YSZ 


-1 


Pore  former 


Fig.  9.  Expansion  surface  response  model  (a)  for  YSZ  and  NiO  particle  size  variation 
with  30wt%  pore-former  and  40wt%  NiO  and  (b)  for  NiO  and  pore-former  ratio 
variation  for  fine  NiO  and  YSZ  powders. 


A  graphical  presentation  to  compare  the  different  relative 
regression  coefficient  estimators  is  given  in  Fig.  10.  100%  corre¬ 
sponds  to  the  value  of  the  constant  or  mean  estimator  (a0). 

Increasing  the  pore-former  concentration  (a2)  as  well  as  the  YSZ 
particle  size  (a4)  decreases  the  expansion  (Fig.  10a).  This  could  be 
related  to  the  porosity  response  effects  (Fig.  lOd),  the  pore-former 
and  the  coarse  YSZ  both  increasing  porosity.  The  coarse  powder 


decreases  the  sinterability,  in  the  opposite  way;  it  increases  the 
porosity  [29].  As  shown  in  Fig.  5,  porosity  and  expansion  after  re¬ 
oxidation  are  correlated.  The  higher  the  porosity,  the  lower  the 
expansion  after  re-oxidation. 

An  interesting  point  is  that  the  NiO  content  in  the  composition 
(ai  in  Fig.  10a)  has  no  influence  on  expansion;  the  value  of  this 
effect  is  due  to  ‘noise’  as  shown  by  the  normal  distribution  of  the 
coefficients  in  Fig.  11a.  From  a3,  coarse  NiO  seems  to  be  the  worst, 
as  it  increases  expansion.  The  pore-former  has  a  positive  quadratic 
influence;  this  is  shown  by  Fig.  9b:  for  fine  YSZ  and  NiO  powders, 
minimal  expansion  does  not  lie  at  the  border  of  the  variables  space 
but  at  around  1 5  wt%  of  pore-former. 

The  wide  range  of  conductivity  values  from  1 6,200  S  cm-1  (com¬ 
position  5)  to  OS  cm-1  (compositions  10  and  16),  induces  large 
estimate  coefficients  of  the  factors  (Fig.  10b).  The  higher  effect 
is  from  the  NiO  content;  it  is  more  than  100%  of  the  constant 
value  (a0).  Further,  conductivity  is  decreased  by  the  pore-former 
content  (a2).  The  conductivity  depends  directly  on  the  volume 
ratio  of  the  metallic  phase  and  the  porosity  of  the  sample  as  seen 
in  the  previous  section.  At  high  NiO  concentration,  the  effect  of 
NiO  particle  size  is  even  stronger  (from  ai3).  Pore-former  con¬ 
tent  together  with  high  NiO  and  coarse  YSZ  particle  sizes  (a23, 
a24)  lead  to  a  very  negative  interaction  effect  on  conductiv¬ 
ity. 

The  surface  quality  decreases  by  adding  pore-former  and 
increasing  the  YSZ  particle  size.  This  is  probably  due  to  the  lower 
sinterability  with  pore-former  and  coarse  YSZ,  to  yield  damaged 
sample  surfaces.  If  both  factors  are  increased,  the  effect  on  surface 
quality  is  even  worse  due  to  the  negative  interaction  estimator 
(a24)  in  Fig.  10c.  The  quadratic  effects  are  important  for  sur¬ 
face  quality  compared  to  the  main  effects:  a  positive  quadratic 
effect  gives  maxima  at  the  borders  of  the  studied  region  for  pore- 
former  and  NiO  particle  size.  For  the  YSZ  particle  size  it  is  at  the 
center  of  the  domain  that  the  maximal  surface  quality  can  be 
found. 

In  summary,  the  major  effect  is  given  by  the  pore-former  addi¬ 
tion,  which  increases  the  as-sintered  porosity  and  reduces  the 
expansion  on  re-oxidation.  The  coarse  YSZ  particles  give  a  similar 
response  due  to  lower  sinterability.  The  NiO  content  only  affects 
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Fig.  10.  Relative  regression  estimator  coefficient  (half-effect)  of  the  quadratic  model  for  (a)  expansion  after  one  RedOx  cycle  at  845  °C,  (b)  conductivity  at  room  temperature 
after  24  h  of  reduction,  (c)  surface  quality  after  sintering  and  (d)  porosity  after  sintering.  All  half-effects  are  divided  by  their  constant  effect  Qq. 


the  conductivity,  with  no  threshold  value  observed.  The  NiO  parti¬ 
cle  size  does  not  show  any  significant  effect  on  the  studied  response 
(the  estimator  coefficient,  03,  is  not  significant),  which  could  be 
due  here  to  the  fact  that  coarse  NiO  was  made  from  agglomerates 
broken  down  during  the  ball-milling  process. 


4.4.  Optimized  compositions 

The  optimal  composition  should  present  a  correct  surface  qual¬ 
ity  (SQresponse)  after  sintering  in  order  to  produce  SOFC  half-cells, 
with  dense  thin  electrolyte  e.g.  for  subsequent  cathode  screen- 


Fig.  11.  Regression  estimator  coefficient  (half-effect)  normal  plot  for  (a)  expansion  after  one  RedOx  cycle  at  845  °C,  (b)  conductivity  at  room  temperature  after  24  h  of 
reduction,  (c)  surface  quality  after  sintering  and  (d)  porosity  after  sintering. 
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Table  7 

Optimized  compositions  fulfilling  the  requirements.  Surface  quality  >  1 .65,  conductivity  at  room  temperature  >  500  S  cm-1  and  expansion  <  0.2%. 


Case 

Compositions 

Response  of  the  models 

NiO  ratio  (wt%) 

Pore-former  (wt%) 

NiO  particle  size 

YSZ  particle  size 

Expansion  (%) 

Conductivity  (S  cirr1 ) 

Surface  quality 

1 

40 

3 

Fine 

Coarse 

0.10 

1554 

2.60 

2 

40 

10 

95%  Coarse 

Coarse 

-0.29 

518 

1.66 

3 

40 

20 

Fine 

Coarse 

-1.48 

837 

1.65 

4 

50 

7.5 

95%  Fine 

Coarse 

0.16 

2402 

2.28 

5 

50 

15 

Fine 

Coarse 

-0.70 

1587 

1.94 

6 

60 

9 

Fine 

Coarse 

0.03 

3247 

2.16 

7 

60 

15 

Fine 

Coarse 

-0.50 

2440 

1.82 

8 

60 

20 

Fine 

Fine 

-0.07 

2807 

2.03 

printing  and  stack  assembling.  It  also  should  possess  sufficient 
conductivity  and  show  an  expansion  after  re-oxidation  lower  than 
0.2%.  The  limitation  of  electrical  conductivity  is  given  by  the  thick¬ 
ness  of  the  anode  support  and  the  distance  between  the  electrical 
contacts  at  the  surface  of  the  anode.  Conductivity  at  800  °C  higher 
than  100  S  cm-1  is  an  adequate  target.  The  variation  of  electrical 
conductivity  of  the  cermet  with  temperature  has  been  measured  by 
Wei  and  Atkinson  [30].  The  electrical  conductivity  at  800  °C  is  about 
20%  of  the  value  at  room  temperature.  Hence  a  limit  of  500  S  cm-1 
measured  on  the  sample  at  room  temperature  can  be  suggested  as 
criterion  (Table  7). 


To  find  an  optimum  composition,  the  lowest  surface  quality  (SQ) 
limit  is  arbitrarily  fixed  at  1 .65,  the  lowest  (room  temperature)  con¬ 
ductivity  value  at  500  S  cm-1  and  the  highest  expansion  limit  after 
re-oxidation  at  0.2%.  The  compositions  that  fulfill  these  conditions 
are  presented  in  four  3D  graphs  with  the  axes  corresponding  to  the 
different  factors.  From  Fig.  12,  it  is  possible  to  find  different  com¬ 
positions  that  are  in  agreement  with  these  requirements.  Negative 
values  of  expansion  are  due  to  the  quadratic  model  fit.  It  is  inter¬ 
esting  to  note  that  appropriate  compositions  can  be  found  over  the 
full  range  of  studied  NiO  content  (40-60  wt%  NiO)  and  YSZ  powder 
size.  The  pore-former  content  can  vary  from  3  to  20wt%  and  the 


Pore  former 


-1  -1 


NiO  proportion 


Fig.  12.  Compositions  where  surface  quality  >  1 .65,  conductivity  >  500  S  cm-1  and  expansion  <  0.2%. 
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NiO  powder  size  from  5  to  100%  fine  fraction.  These  compositions 
will  in  future  be  tested  as  anode  supports  for  SOFC. 


5.  Conclusion 

A  design  of  experiment  approach  was  efficiently  applied  to  the 
microstructure  optimization  of  the  Ni/YSZ  anode  support  for  RedOx 
stability  enhancement.  Quadratic  models  give  a  significantly  better 
experimental  fit  than  linear  models. 

The  major  advantage  of  this  work  compared  to  previous 
studies  is  the  wide  domain  of  experimental  space  (81  possible 
compositional  combinations)  covered  with  a  reduced  number  of 
experiments  (25  variance-optimized  compositions  and  46  tested 
samples),  with  a  method  giving  the  statistics  to  extract  individual 
and  combined  effects  of  the  different  experimental  variables  on  the 
study’s  final  target,  i.e.  the  limited  re-oxidation  strain  of  the  Ni-YSZ 
anode  support. 

The  RedOx-safe  expansion  limit  lower  than  0.2%  is  reached  only 
for  samples  with  as-sintered  porosity  at  least  higher  than  35%,  and 
rather  around  50%,  for  the  parameter  space  covered  in  this  study. 
The  major  factors  lowering  the  expansion  are  correlated  to  the  ones 
increasing  porosity:  pore-former  fraction  and  YSZ  powder  size. 
Coarse  zirconia  powder  has  lower  sinterability  and  so  increases 
final  sample  porosity. 

Conductivity  is  influenced  principally  by  the  NiO  content  and 
secondly  by  the  as-sintered  porosity. 

As  a  result  of  the  approach,  eight  compositions  reaching  expan¬ 
sion  smaller  than  0.2%,  conductivity  higher  than  500  S  cm-1  and 
good  sample  surface  quality  are  defined.  Forthcoming  work  will 
consist  in  selecting  several  of  these  compositions  to  be  produced 
by  tape  casting  and  electrochemically  tested. 


Acknowledgments 


Table  B1 

Density,  molar  mass  and  molar  volume  of  the  components  of  the  studied  composite. 


Ni 

NiO 

8YSZ 

Pl  (gem-3) 

8.89 

6.67 

5.90 

M/(gmol-1) 

58.71 

74.71 

127.87 

V/( cm3  mol-1) 

6.604 

11.20 

21.67 

Weight  to  molar: 

wM2 


wM2  +  (1  -  w)Mi 
Molar  to  weight: 

_  nMi 

W=  nM,  +  (1  -  n)M2 

Weight  to  volume: 

v=  Wp2 

wp2  +  (1  -  w)pi 

Volume  to  weight: 

vpi 

w  = - — - 

VP\  +(1  -  V)p2 

Volume  to  molar: 

vV2 

n~  w2  +  (1  —  t/)V, 
Molar  to  volume: 
nVi 


nVi  +(1  -n)V2 

Volume  oxidized  to  volume  reduced: 

_  _ v 0X^1, red _ 

fox^l,red  +  (1  —  Vox)V\  ,ox 

Weight  oxidized  to  weight  reduced: 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
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Appendix  A.  Supplementary  data 


wred 


Wox^i  ,red 

WoxMired  +  (1  —  Wox)Ml,ox 


(B9) 


Including  porosity  fraction  in  oxidized  state  (same  for  the  reduced 
state): 


Pi, ox  +P2,ox  +Ppore,ox  —  1 


Pi,  ox 
fl,ox 


P2,ox 
^2  ,ox 


P pore,  ox 


Porosity  after  reduction  (for  a  constant  sample  volume): 


(B10) 


Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  doi:  10.1 01 6/j.jpowsour.201 0.07.092. 

Appendix  B 

Here  the  derivation  is  given  to  calculate  the  different  fractions 
(weight,  volume  and  molar)  in  a  composite  (w:  weight  fraction,  v: 
volume  fraction  in  the  solid  phase,  n:  mole  fraction,  p,- :  volume  frac¬ 
tion  including  porosity)  [2,31].  The  composite  is  based  on  phase  1 
and  2.  When  the  subscript  is  not  written,  the  fraction  is  for  compo¬ 
nent  1 .  In  the  case  of  Ni/NiO-YSZ  composite  the  subscripts  are  1 ,  red 
for  Ni,  l,ox  for  NiO,  2  for  YSZ,  and  pore  for  porosity.  The  density  is 
p,  the  molar  volume,  V  and  the  molar  mass,  M.  The  sample  volume 
is  considered  as  constant  during  the  state  change  (from  oxidized  to 
reduced  state). 

Density  of  the  composite: 

PiP2 

Pth  wp2  +  (l-w)pi 


P pore, red  —  ~ AVrec/V l,ox(l  ~  P pore, ox)  +  P pore, ox  (B11) 

with  AVred  =  -  1  (B12) 

v1,ox 
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